Unconventional Vortex Dynamics in Mesoscopic Superconducting Corbino Disks 
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The discrete shell structure of vortex matter strongly influences the flux dynamics in mesoscopic 
superconducting Corbino disks. While the dynamical behavior is well understood in large and in 
very small disks, in the intermediate-size regime it occurs to be much more complex and unusual, 
due to (in)commensurability between the vortex shells. We demonstrate unconventional vortex 
dynamics (inversion of shell velocities with respect to the gradient driving force) and angular melt- 
ing (propagating from the boundary where the shear stress is minimum, towards the center) in 
mesoscopic Corbino disks. 
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PACS numbers: 74.25.Qt, 74.78.Na, 68.35. Af 

Investigation of mesoscopic and nano-structured su- 
perconductors has become one of the major topics of 
experimental and theoretical studies on vortex matter 
in superconductors over the last two decades. The 
unique possibility to investigate the interplay between 
the vortex-vortex and vortex-boundary interactions, the 
size-, shape-, and symmetry- induced effects makes meso- 
scopic superconductors very attractive from both the fun- 
damental research and possible fluxonics applications, 
i.e., manipulating single flux quanta similar to electrons 
in nano- and microelectronics. Thus it was demonstrated 
theoretically P, 0, II H and experimentally @, i, 0, Q 
that in mesoscopic superconducting disks shape effect 
leads to the formation of circular-symmetric giant vortex 
states d, @ or concentric shells of vortices @, [3 ■ Further- 
more, symmetry-induced vortex-antivortex "molecules" 
were theoretically predicted [10] in mesoscopic squares 
and triangles. On the other hand, the dynamics of vor- 
tex matter in disks can be studied using a Corbino setup. 
In the Corbino geometry an applied current produces a 
gradient in the Lorentz force, i.e., introduces a shear driv- 
ing force between the rings (in large disks [ll|) or shells 
(in mesoscopic regime [12| ) of vortices. This provides 
the unique opportunity to study various dynamical ef- 
fects related to vortex motion, e.g., the onset of plastic- 
ity, channeling [13], vortex friction [l4|], etc. Moreover, 
the dynamics of self-organized vortex matter in meso- 
scopic disks has many common features to, e.g., atomic 
matter, charged particles in Coulomb crystals, vortices 
in rotating Bosc-Einstein condensates, magnetic colloids, 
synthetic nanocrystals, etc. [IH fltl [l7| . o r charged balls 
diffusing in macroscopic Wigner rings [181 ]. 

Several early studies on vortex matter in Corbino disks 
focused on the transition from elastic to plastic vortex 
motion [ll|, 0, [2(|. Indeed, the radial current density 
j in a Corbino disk decays as Io/r along the radius re- 
sulting in a stronger Lorentz force Fl driving vortices 
near the disk center than near the edge. For small driv- 
ing currents, the local shear stress is small and the whole 
vortex pattern, although elastically deformed, moves as a 
rigid body. For applied currents larger than some critical 
value I c a strong spatially inhomogeneous stress breaks 



up the vortex solid, and concentric annular regions (in 
large disks) or shells (in mesoscopic disks) move with dif- 
ferent angular velocities. The different dynamical phases 
(e.g., elastic motion, shear-induced plastic slip) of vor- 
tex motion can be detected by measuring voltage pro- 
files [ll[. In large Corbino disks, the onset of plastic- 
ity was theoretically analyzed within a continuous model 
20] and using molecular dynamics (MD) simulations of 
interacting vortices at T = [2l|. The magnetization 
and different scenarios of magnetic-flux penetration in a 
thin Corbino disk were studied [22j using the critical- 
state model. The dynamics of vortex shells in small 
mesoscopic Corbino disks containing up to 40 vortices 
(cf. Ref. Q) arranged in two- and three-shell configura- 
tions, was studied using MD simulations for T — and 
T > 12]. It was shown that the critical current I c had a 
non-monotonic dependence on the magnetic field h. The 
dynamical instabilities associated with a jump in I c (h) 
were revealed. The found unusual behavior was shown to 
be related to a "structural transition", i.e., an inter-shell 
vortex transition with increasing h. In the presence of 
pinning, thermally-activated externally-driven flux mo- 
tion was analyzed for T > 0. The investigation of the an- 
gular thermal melting explained the experimentally ob- 
served two-step melting transition in Corbino disks [111 ]. 

Here, we study the dynamics of vortex shells in larger 
mesoscopic Corbino disks containing several shells. The 
essential qualitative difference to the earlier studied case 
fl2| is the following: While in small disks vortex shells 
are nearly perfectly circular, due to the strong confine- 
ment, in larger disks the vortex-vortex interaction be- 
comes more important. The dynamics of vortex shells is 
influenced, to a great extent, by the interplay between 
the gradient Lorentz force and the (in)commensurability 
effects between the numbers of vortices in adjacent (and 
even more remote!) vortex shells. Note that all these 
factors are sensitive to the radius of the disk. As a re- 
sult of the complex interplay of the above five factors, 
i.e., a) confinement, b) the inter- vortex interaction, c) 
the gradient Lorentz force, d) (in)commensurability, and 
e) the size of the disk, the system of vortex shells in 
large Corbino disks demonstrates very rich and unusual 
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dynamical behavior. Thus we found: (i) unconventional 
angular melting, i.e., when the melting first occurs in re- 
gions where the shear stress is minimum, and (ii) uncon- 
ventional shell dynamics when a shell that experiences a 
weaker Lorentz force moves faster than the adjacent shell 
driven by a stronger force. We expect that our results can 
be useful for understanding and predicting the dynamics 
of complex systems of interacting particles confined to 
an external potential such as colloids, vortices in Bose- 
Einstein condensates, charged particles, etc., or even for 
better understanding of mechanisms of motion of some 
biological objects. 

Theory and simulation. — We place a Corbino disk 
which has thickness d and radius R in a perpendicu- 
lar external magnetic field Hq. An external current 
flows radially from the center to the edge of the disk 
and results in the inhomogeneous sheath current density 
j(p) — I ex t/2npd. The Lorentz force (per unit length) 
acting on vortex i, $oj X z, resulting from the external 
current is: 

f * -^d-—- foI °—' (1) 

where pi is the vortex position, Ti — pi/R, and z is 
the unit vector along the magnetic field direction. Here 
the unit of force f = ^1/2tt^ RX 2 = 4n n £ 2 H 2 / R and 
= P-o-h-Iext/^o are expressed in terms of the coher- 
ence length £, the magnetic field penetration depth A 
(A = X 2 /d is the effective penetration depth in a super- 
conductor of thickness d), and the thermodynamic H c 
and the upper H C 2 critical magnetic fields. In a thin 
superconducting disk such that d < £ <C R <C A, the 
vortex-vortex interaction force /™ and the force of the 
vortex interaction with the shielding currents and with 
the edge // can be modelled respectively by [H, [24], [2f| 

rr = fot - ^ ,f*~% ) , (2) 

i,k \\ r i-r k \ \rln-r k \ J 

# = /o(i^f -h)n, (3) 

where h = (Hq/2H C 2)(R/^) 2 is the dimensionless ap- 
plied magnetic field, and L is the number of vortices, 
or the vorticity. Our numerical approach is based on 
the Langevin dynamics algorithm, where the time in- 
tegration of the equations of motion is performed in 
the presence of a random thermal force. The over- 
damped equation of motion (see, e.g., [2(1 [2?|]) becomes: 
VVi = fi = fr+fr+fi+fi+fi- Here ft is the driv- 
ing force (Eq. ((!])), f? p is the force due to vortex-pin in- 
teraction, and fj is the thermal stochastic force, obeying 
the condition: (f a ,i(t)fp,j(t')) = 2r)5 al3 5 i] S(t- t')k B T, 
where Greek and italic indices refer to vector components 
and vortex labels, r] is the viscosity which is set to unity. 
The ground state of the system is obtained by minimiza- 
tion of the London free en erg y, taking into account the 
vortex cores contribution [281 ]. with respect to: (i) dif- 
ferent L, and (ii) different vortex configurations for the 




0.000 0.002 0.004 0.006 



'o 

FIG. 1: (Color online) The angular velocities uji of different 
vortex shells as a function of the external driving current 7o. 
The critical values 7 c i2 and I C 23. correspond to the unlock 
transition between the shells S 1 -S 2 and S 2 -S 3 . The two-step 
angular melting starts at the center of the disk where the 
Lorentz force is maximum. The inset: The shell configuration 
(1,5,12,18) for the vorticity L = 36 and h = 45. The dashed 
grey circles showing different shells are guides for eye. 

same L (see Refs. @, H5) by performing the stimulated 
annealing simulation (SAS) that simulates field-cooling 
experiments. 

Angular melting: Commensur ability vs. the Lorentz 
force gradient. — In a Corbino disk the gradient of the 
Lorentz force is maximum near the center, and thus one 
can expect that with increasing applied current the grow- 
ing shear stress first breaks the vortex solid near the cen- 
ter and then, with further increasing the applied current, 
the melting process [29| propagates from the center to 
the peripheral regions. Indeed, this scenario of angular 
melting of the vortex solid in Corbino disks was observed 
in the experiment [ill ] and analyzed theoretically (e.g., 
0, H3, El| ) ■ Similarly, in mesoscopic disks vortex shells 
first unlock near the center with increasing either applied 
current or temperature [l2l ] . Fig. 1 shows the angular ve- 
locity of different vortex shells u>i for the configuration 
(1,5,12,18) as a function of the applied current Iq. At 
low Iq, the whole vortex configuration rotates as a rigid 
body while at some critical value I c i2 the first inner shell 
splits off and starts to rotate separately with a higher 
angular velocity uj\. At a higher driving current J C 23, the 
Lorentz force gradient becomes sufficient for unlocking 
the second shell. For vortex configuration with a larger 
number of shells, the angular melting thus occurs gradu- 
ally through the consequent unlocking of the shells from 
the center towards the periphery. 

It is easy to understand that the above "classical" be- 
havior is the only possible scenario for angular melting of 
the vortex solid in large (macroscopic) Corbino disks. In- 
deed, in large disks (i.e., where the intervortex distance 
a <C R) vortices form a uniform hexagonal Abrikosov 
lattice. It is clear that the shear modulus of such lat- 
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FIG. 2: (Color online) Unconventional angular melting of 
vortex shells for the configuration (1,6,12,19): the two-step 
melting starts near the boundary where the Lorentz force is 
minimum. The inset shows the shell configuration (1,6,12,19) 
for the vorticity L = 38 and h = 50. 

tice is the same everywhere (except for the regions in 
the vicinity to the boundary). Thus, the strongest shear 
stress which occurs near the center of the disk breaks 
up the lattice first. However, in case of vortex shells in 
mesoscopic disks the situation turns out to be very differ- 
ent: the shear modulus becomes a function of the radius 
due to (in)commensurate number of vortices in differ- 
ent adjacent shells. (Moreover, the shear modulus can 
also depend on the angle if the circular symmetry of the 
vortex-shell configuration is broken by defects.) Thus we 
can expect that the angular melting in a system of vortex 
shells will happen in a different way if the shear modulus 
has a strong radial dependence, i.e., stronger than the 
difference in Lorentz force between adjacent shells. 

Because of the complexity of the system it is not possi- 
ble to find the shear modulus analytically. However, we 
know that the shear modulus will be much smaller for 
incommensurate shells than for the commensurate case 
where "magic number configurations" [30 . l3l| make the 
vortex system more rigid. If the inner shells are incom- 
mensurate, the Lorentz force gradient will facilitate the 
"classical" scenario of angular melting of the vortex solid 
and the unlocking of shells will occur at an even lower 
value of the driving current Iq. In contrast, if the outer 
shells are incommensurate, the friction between them can 
be so small that they first unlock, leading to unconven- 
tional angular melting, i.e., from the outer edge of the 
vortex solid inwards. 

We analyzed different ground-state and metastable 
vortex-shell configurations for different number of vor- 
tices L and found different scenarios of angular melting 
in a mesoscopic Corbino disk. The most unusual way 
of melting is when it starts from the outermost shell 
and propagates towards the center of the disk, i.e., op- 
posite to the "classical" way. An example of such un- 
conventional melting is presented in Fig. 2. Note that 
it is realized for the ground-state energy configuration 
(1,6,12,19) for L = 38. The circular-like incommensurate 
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FIG. 3: (Color online) Ui(Io) for the configuration 
(1,6,12,19,25,30). The system displays unconventional vortex 
dynamics in the region 7 C 234 ~ 0.004 < 7o < J C 34-43 ~ 0.0095 
where the shell S 4 (with 25 vortices), that experiences a 
weaker Lorentz force, rotates faster than the adjacent shell 
5 3 (19 vortices) driven by a larger Lorentz force. The inset 
shows the shell configuration (1,6,12,19,25,30) for the vorticity 
L — 93 and h — 109. The S 3 -shell consists of a hexagon-like 
part (CP) commensurate with the the hexagonal core (1,6,12) 
and a "circular-like" part (ICP) (with a defect) incommensu- 
rate with the core. 

shell with 19 vortices easily slides with respect to the 
adjacent shell with 12 vortices even for very low driv- 
ing currents I > I C 23 while the hard "core" (1,6,12) can 
sustain much stronger shear stress and the inner shell 
becomes unlocked at I c i2 ^> I C 23 (Fig. 2). 

Unconventional dynamics of vortex shells. — The in- 
terplay between (in)commensurability effects of vortex 
shells, on the one hand, and the inhomogeneous Lorentz 
force, on the other hand, can lead to even more strik- 
ing and unexpected behavior. Thus, we demonstrate 
here that these effects can result in a very unusual ( "un- 
conventional" ) dynamics of vortex shells in a mesoscopic 
Corbino disk. In particular, we found that in the region 
of the multistep transition from a "rigid body" rotation 
to individual rotation of vortex shells, an inversion of 
the angular velocities with respect to the gradient of the 
Lorentz force can occur. For unlocked shells we expect 
that, due to the \ jr dependence of the driving force, the 
angular velocity of the different shells is a monotonically 
decreasing function of the radius. This is clearly seen 
in Figs. 1 and 2. Very surprisingly, we found situations 
were this rule is broken as demonstrated in Fig. 3 for 
the L = 93 configuration (1,6,12,19,25,30). The incom- 
mensurate shell S with 25 vortices rotates faster than 
the adjacent shell S 3 (which is closer to the center than 
S 4 ) for 7 c2 34 ~ 0.004 < I < I c3 4-43 ~ 0.0095. This 
puzzling behavior is a consequence of the interplay be- 
tween the Lorentz force, (in)commensurability effects and 
shell defects. For the above range of Iq, the vortex-shell 
configuration is only partially split. The strong commen- 
surate inner-shell "core" rotates as a rigid body while 
incommensurate outer shells seem to rotate separately. 
However, their motion turns out to be correlated in a 
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rather complex manner. Namely, the S^-shell with 19 
vortices has a "defect", i.e., one extra vortex (indicated 
by a blue circle in the inset of Fig. 3) which makes this 
shell incommensurate with the hexagon-like core (1,6,12). 
However, due to the inter-shell vortex-vortex interaction 
a part of this shell squeezes and adjusts itself to the core 
thus making the other part of the shell (i.e., incommensu- 
rate, containing the defect) stretched. The commensurate 
part (CP) of the S^-shell (the hexagon-like part of shell 
S 3 - see inset of Fig. 3) is locked to the core while the 
incommensurate part (ICP) (the "circular" part of shell 
S 3 - see inset of Fig. 3) is unlocked and slides with re- 
spect to the core. Thus the motion of the 5' 3 -shell with 
respect to the core occurs in the form of a compression- 
decompression wave propagating in the direction opposite 
to the direction of the rotation. This backward wave mo- 
tion reduces the average angular velocity of the 5 3 -shell 
by a factor of 18n/19 where n — 2irrs3/Xcp-icp an d 
^cp-icp is the wavelength of the CP-ICP compression- 
decompression wave [Hj]. At the same time, the incom- 
mensurate S^-shell with 25 vortices, due to the friction 
with the CP of the S^-shell, is accelarated by the faster 
rotation of the core. Note that this counter-intuitive be- 
havior is only possible when a very delicate balance be- 
tween the vortex- vortex interaction and the Lorentz force 
is established. When the latter becomes large enough 
(i.e., for Io > 0.01) the usual vortex-shell velocity dis- 
tribution (i.e., monotonic decreasing) is recovered (see 



Fig. 3). Note that at Iq w 0.01 the angular velocity of 
the core also has a feature (kink) indicating the complete 
unlocking between the shells S 3 , S 4 and the core (which 
itself breaks apart only for very large values of driving, 
Io > Ici2 ~ 0.0325). 

Conclusions. — We found that the interplay between 
the (in)commensurability effects, the vortex-vortex in- 
teraction and the gradient Lorentz force in a mesoscopic 
Corbino disk can result in a very unusual and counter- 
intuitive dynamics of vortex shells. Thus we demon- 
strated: (i) unconventional angular melting when, with 
increasing driving current, the melting first occurs near 
the boundary where the shear stress is minimum and 
propagates towards the center, in an opposite way to the 
"classical" scenario of conventional angular melting in a 
Corbino disk; (ii) unconventional shell dynamics when a 
shell that experiences a weaker Lorentz force moves faster 
than the adjacent shell driven by a stronger Lorentz force. 
These examples of strongly non-linear dynamical behav- 
ior can be useful for understanding and prediction of dy- 
namics of other complex interacting systems. 
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